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Abstract We recently presented a model for the cluster formation efficiency (CFE), 
i.e. the fraction of star formation occurring in bound stellar clusters. It utilizes the 
idea that the formation of stars and stellar clusters occurs across a continuous spec- 
trum of ISM densities. Bound stellar clusters naturally arise from the high-density 
end of this density spectrum. Due to short free-fall times, these high-density regions 
can achieve high star formation efficiencies (SFEs) and can be unaffected by gas ex- 
pulsion. Lower-density regions remain gas-rich and substructured, and are unbound 
upon gas expulsion. The model enables the CFE to be calculated using galactic- 
scale observables. I present a brief summary of the model physics, assumptions and 
caveats, and show that it agrees well with observations. Fortran and IDL routines for 
calculating the CFE are publicly available at http://www.mpa-garching.mpg.de/cfe 



1 The clustered nature of star formation 

Studies of star formation in the solar neighbourhood, in the Milky Way as a whole, 
and in external galaxies, have shown that some fraction of star formation occurs 
in unbound associations, while the remainder results in bound stellar clusters (e.g. 
[UEl). The question thus arises which physical mechanisms drive the star formation 
process to result in either outcome. How do bound stellar clusters form, and what 
fraction of all star formation do they represent? Observations suggest a range of 
answers to the latter question - it seems that the fraction of star formation occurring 
in bound stellar clusters changes with the galactic environment [3] 0J [5] |6) . This 
provides a potentially useful clue to the origins of stellar clustering, and presents a 
challenge to the classical view in which all stars form in clusters but only a minor 
fraction remains bound after gas expulsion (e.g. J7)). 
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Another important insight is that star formation occurs over a broad and continu- 
ous range of densities O, without any obvious separation between star formation in 
unbound associations or bound clusters. Instead, protostellar cores and young stellar 
objects follow the hierarchical structure of the interstellar medium (ISM) (9J. It has 
recently been shown that gas-poor and virialised stellar structure may arise naturally 
at the high-density end of this density spectrum ifTOlfTTIl . This may provide a nat- 
ural explanation to the observation that cluster formation is inefficient lfl2l : if only 
some fraction of the star-forming regions manages to collapse and form bound sys- 
tems, that same fraction will emerge as unembedded, virialized and bound clusters 
without any sign of expansion due to gas expulsion (e.g. |[T3l ). 

In 03], we recently presented a new theoretical framework for the formation of 
bound stellar clusters. By integrating the star formation efficiencies (SFEs) and local 
bound fractions of star-forming regions over the density spectrum of the ISM, our 
model is used to quantify the fraction of all star formation that occurs in bound stel- 
lar clusters (i.e. the cluster formation efficiency or CFE). Fortran and IDL routines 
for calculating the CFE as a function of galaxy properties are publicly available at 
http://www.mpa-garching.mpg.de/cfe [14|. In these proceedings, I summarize the 
model and discuss the underlying assumptions and caveats. 



2 The fraction of star formation occurring in bound clusters 

The theory of the CFE is derived and applied in detail in lfl4l . In summary, it covers 
the following physical mechanisms and underlying assumptions. 

1. The starting point of the model is the overdensity probability distribution func- 
tion (PDF) of the ISM. This PDF is assumed to follow a log-normal with median 
and dispersion set by the Mach number, and describes the distribution of density 
contrasts with respect to the mean density in a turbulent ISM (see e.g. lfT31[T6l ). 

2. The overdensity PDF is written as a function of galaxy properties (gas surface 
density E g , angular velocity £2, and Toomre Q 021) by assuming that star for- 
mation occurs in a gas disc that obeys hydrostatic equilibrium (see below). 

3. At each density the local SFE is calculated by assuming that the fraction of 
the gas that is converted into stars per free-fall time is approximately constant 
lfT8l[T31 . Depending on the density, star formation continues until (1) the gas is 
exhausted, (2) pressure equilibrium is reached between the turbulent ISM and 
(supernova and/or radiative) feedback, or (3) the moment of evaluating the CFE. 

4. The local SFE is related to the local fraction of stars that remains bound upon in- 
stantaneous gas expulsion using a numerical simulation turbulent fragmentation 
|fl9l . Assuming that protostellar outflows do not unbind bound stellar clusters, 
this provides the naturally bound fraction of star formation at each density. 



5. Following the Spitzer formalism for tidal shocks 1201 . the model includes tidal 
perturbations by density peaks in the star-forming environment (the cruel cradle 
effect, see 12T1[TT1 ). which destroy stellar structure below a certain environmen- 
tally dependent, critical density. 
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6. The CFE is obtained by integrating the naturally bound fraction of star formation 
over the density range of the PDF where structure survives the cruel cradle effect 
(reflecting bound cluster formation), and dividing it by the integral of the SFE 
over the entire density range of the PDF (reflecting all star formation). 

Given these model components, the following caveats should be kept in mind. 

1 . Magnetic fields are only included to first order by using the magnetic-to-thermal 
pressure ratio, which is specified with an optional model parameter. This changes 
the dispersion of the log-normal overdensity PDF [16]. Note that the default form 
of this PDF is already consistent with weak magnetic fields. 

2. The assumption that star formation occurs in a gas disc that obeys hydrostatic 
equilibrium may not be consistent with starburst galaxies. However, the energy 
dissipation that is required to cool gas and form stars also drives the formation 
of a disc, and hence the spatial distribution of star-forming regions in a starburst 
should be expected to follow a disc-like morphology l22l . 

3. Star formation in intermediate-density regions is halted by supernova (or alter- 
natively radiative) feedback, which may not be appropriate. The efficiency of 
different feedback mechanisms has been extensively discussed in the literature, 
and likely varies with spatial scale or density. However, the description of feed- 
back in this model satisfies its purpose of truncating star formation on a timescale 
that is broadly consistent with observations (TJ. 

The above model can be used to calculate the CFE as a function of {E g ,Q,Q}, 
which is reduced to a single-parameter problem by assuming a single value of 
Q = 1.5 and relating Q to E g as in lfT31 . The resulting relation between the CFE 
(or r, 11231 ) and the star formation rate density Zsfr x ll24l is shown in Fia- 
ureQ] together with compiled observations from the recent literature. The agreement 
between theory and observations is remarkable, especially considering that the typ- 
ical error margins on the observations are ~ 0.3 dex, and noting that an additional 
uncertainty of 0.3-0.5 dex is introduced by using the relation Xsfr x £g 4 to convert 
the model r(E g ) to .T(Zsfr)- A good fit to the model relation is given by 



where Zsfr.o = ^SFr/M© yr~ 1 kpc~ 2 is the star formation rate density of the galaxy. 
This fit is for one particular, 'typical' parameter set and should therefore only be 
used for rough estimates. Because it assumes the power law form of the Schmidt- 
Kennicutt relation £sfr x -Eg 4 ll24l . any scatter around that relation is carried over. 
See lfl4l and our publicly available routines for a more detailed modelling. 

The good agreement between model and observations warrants further testing 
using Gaia and ALMA (see |[T4| for an extensive discussion of the possibilities). 
However, our theoretical framework also contains several components that require 
to be constrained further. While a global theoretical picture of cluster formation 
seems to be emerging, the details of stellar clustering remain to be understood. 
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Fig. 1 CFE as a function of 
•£sfr- Symbols denote ob- 
served galaxies with 1 o~ error 
bars and indicate the samples 
from [3] (blue diamonds), 
|4| (green triangles), and |5| 
(red squares). The black cross 
indicates the integrated CFE 
of all dwarf galaxies from the 
sample of (6), with a surface 
density range indicated by 
the horizontal error bar. The 
solid curve represents the 
modelled relation for typical 
disc galaxies, but can vary 
for different galaxy and ISM 
properties. The dashed curve 
shows the fit of equation {T}, 
and the dotted line represents 

the original fit by (3)- ^sfr L m o J t ' k P c 
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